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Abstract. We have analysed the ISO-SWS spectrum of 
Jupiter in the 12-16 /im range, where several hydrocar- 
bons exhibit rovibrational bands. Using temperature in- 
formation from the methane and hydrogen emissions, we 
derive the mixing ratios (q) of acetylene and ethane at 
two independent pressure levels. For acetylene, we find 
q = (8.9±J;J) x 10" 7 at 0.3 mbar and q = (l.lio.?) x 10~ 7 
at 4 mbar, giving a slope — din q/ din P = 0.8 ±0.1, 
while for ethane q = (1.0 ± 0.2) x 10~ 5 at 1 mbar and 
q = (2.6lg g) x 10" 6 at 10 mbar, giving -din q/ din P = 
0.6 ± 0.2. The ethane slop e is co nsistent with the predic- 
tions of Gladstone et al. ( 1996 ), but that predicted for 
acetylene is larger than we observe. This disagreement is 
best explained by an overestimation of the acetylene pro- 
duction rate compared to that of ethane in the Gladstone 
et al. ( [1996D model. 

At 15.8 /mi, methylacetylene is detected for the first 
time at low jovian latitudes, and a stratospheric column 
^ ' density of (1.5 ± 0.4) x 10 15 molecule cm" 2 is inferred. 



We also derive an upper limit for the diacetylene column 
density of 7 x 10 13 molecule cm" 2 . 

Key words: Planets and satellites: Jupiter -Infrared: so- 
lar system 



1. Introduction 

Hydrocarbons in Jupiter are produced in a series of 
chemical pathways initiated by the photolysis of methane 
in the upper stratosphere. Vertical transport, mainly tur- 
bulent diffusion, redistributes the molecules throughout 
the stratosphere and down to the troposphere, where they 
are eventually destroyed. Hydrocarbons, and particularly 
the most stable of them, are therefore good tracers of the 
upper atmospheric dynamics. In addition, as they strongly 
contribute to the atmospheric opacity in the UV and IR, 
hydrocarbons act as major sources and sinks of heat, 
thereby participating to the stratospheric dynamics. 



All these reasons have strongly motivated theoretical 
studies of the jovian stratospheric photochemistry (Stro- 
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bel |1969t Yung fc Strobel |1980| ; and most recently Glad- 
stone et al. 1996 ). Although nowadays very detailed, these 
models still need to be constrained by observations of mi- 
nor species. However, prior to the ISO mission, only two 
molecules, acetylene (C2H2) and ethane (C2H6), had been 
detected, except in the auroral zones where several other 
minor speci es (C 2H4, C3H4 and C§Kq) have been observed 
(Kim et al. 1985 ). Mean stratospheric mole fractions have 
been inferred for C2H2 and C2H6 by various authors, but 
no precise information on their vertical variations was 
made available. In this paper, we analyse the ISO-SWS 
spectrum of Jupiter between 7 and 17 /zm, in order to de- 
termine the vertical distributions of C2H2 and C2Hg, and 
to search for more complex (C3 and C4) molecules. Sect. 2 
presents the observations. Our analysis of the spectrum is 
presented in Sect. 3. The results are compared with pre- 
vious observations and theoretical predictions in Sect. 4. 

2. Observations 

Descriptions of the Infrared Space Observatory (ISO) 
and of the Short Wavelength Spectrometer (SWS) can be 
found respectively in Kessler et al. (1996) and de Graauw 
et al. (1996). A preliminary analysis of the Jupiter SWS 
spectrum can be found in Encrenaz et al. (1996). New 
ISO-SWS grating observations of Jupiter were obtained 
on May 23, 1997 UT using the AOT 01 observing mode. 
These observations have an average spectral resolution of 
1500, and range from 2.4 to 45 /zm. However, the useful 
range is limited to 2.4-17 /im, due to partial saturation 
at longer wavelengths. The instrument aperture, 14 x 20 
arcsec 2 at A < 12.5 /mi and 14x27arcsec 2 at A > 12.5 /tm, 
was centered on the planet with the long axis aligned per- 
pendicular to the ecliptic, thus roughly parallel to the N-S 
polar axis. It covered latitudes between 30°S and 30°N, 
and ±20° longitude range from the central meridian. The 
absolute flux accuracy is ~ 20%. 

Instrumental fringing generates a spurious signal be- 
tween 12.5 and 17 /mi, which amounts to ~10% of the con- 
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Fig. 1. Comparison between ISO-SWS spectra (solid line) 
and synthetic spectra (dotted line) in the CH4 band 
(upper panel) and the H2 S(l) line (lower panel) 
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Fig. 2. Vertical profiles of C2H2 (upper panel) and C2H6 
(lower panel) used for the calculation of the synthetic spec- 
tra of Fig. H and Fig. [l| The best-fit profiles are shown as 
s olid lines. The triangles are from the Gladstone et al. 
( |1996| ) model. 



tinuum level. This parasitic signal was for the most part 
removed by fitting the detector relative response function 
to the observed spectrum and then dividing it out. Resid- 
ual fringes were further removed by selective frequency 
filtering. 



3. Analysis 

We analysed the ISO-SWS spectrum using a standard 
line-by-line radiative transfer code adapted to Jupiter's 
conditions. We included the molecular absorptions by 
NH 3 , CH 4 , C 2 H 2 , C 2 H 6 , and CH 3 C 2 H using the spectro- 
scopic param eters of the GEISA97 databank (Jacquinct- 
Husson et al. 



1999). We also considered C4H2 absorption, 



using a linelist provided by E. Arie (private communica- 
tion) and band intensities from Koops et al. ( 1984 ). Spec- 
troscopic parameters for the H 2 S(l) line were calcu lated 
using molecular constants from Jennings et al. (1987) and 
Router & Sirota ( 1994 ). The H 2 -He collision- induced con- 
tinuum was calculated following the work of Borysow et 
al. ( |1985| , |1988D - Th e NH 3 vertical distribution was taken 
from Fouchet et al. (|2000[). 



3.1. Temperature profile 

We first calculated synthetic spectra in the region of 
the CH4 1/4 band. For CH4, we used a deep mixing ratio of 
2.1 x 10~ 3 (Niemann et al. 1998| ) and the vertical profile de- 
rived by Drossart et al. (1999) from the CH 4 fluorescence 
emission at 3.3 /Ltm. The band allows one to retrieve 4 
independent points on the temperature profile between 35 
and 1 mbar. We also generated synthetic spectra of the H 2 
S(l) rotational line at 17 \xm. This line probes a broad at- 
mospheric layer at 3-30 mbar. The ortho-to-para ratio of 
H 2 was assumed to follow local thermodynamical equilib- 
rium. The stratospheric temperature profile was adjusted 
in order to best match the absolute emission in the CH 4 
band, and the line-to-continuum ratio of the H 2 S(l) line. 

In practice, starting with the temperature profile mea- 
sured in-situ by the Galileo Probe (Seiff et al. 1998[ ), it was 
necessary to cool it by 2 K between 30 and 5 mbar. At pres- 
sures lower than 5 mbar, the initial profile was warmed by 
4 K up to 165 K. The 2/4 band is also somewhat sensitive 
to the temperature around the 10-^ibar pressure level. We 
found that the temperature remains constant within a few 
degrees between 1 mbar and 1 /zbar, as already observed 
by Seiff et al. (1998). At pressures lower than 1 /ibar, we 
adopted the measurements of Seiff et al. ( 1998 ), vertically 



smoothed in order to remove oscillations due to gravity 
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Fig. 3. Comparison near 13.7 /mi between the ISO-SWS 
spectrum (upper panel) and three synthetic spectra (lower 
panel) calculated with the C2H2 vertical distributions of 
Fig. Solid line gives the best fit. 

waves, noting that our measurements are essentially in- 
sensitive to this pressure range. The fit to the CH4 and 
H2 emissions is presented in Fig. [Ij. 

The 20% uncertainty on the absolute flux calibration 
directly results in an uncertainty of ±2K on the temper- 
ature profile inferred from the CH4 emission. This un- 
certainty partly explains the minor disagreement between 
the modelled and observed H2 line. Indeed, our model, 
while giving an optimum fit to the CH4 emission, slightly 
(5-10%) overpredicts the observed linc-to-continuum ratio 
of the S(l) line. We also note that the H 2 ortho-to-para 
ratio could differ from the thermal equilibrium value, as 
observed in the troposphere by Conrath et al. ( 199§| ). For 
example, a synthetic spectrum calculated with a constant 
para fraction of 0.34, corresponding to the thermal equi- 
librium value at 115 K, would fully reconcile the CH4 and 
H2 measurements. 
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Fig. 4. Same as Fig. |for C 2 H 6 at 12 

obtained with three distinct C2H2 vertical distributions 
(Fig. ||). All three profiles reproduce the emissions due 
to the P-,Q-, and R-branches of the main 1/5 C2H2 band 
(Fig. ||). However, only one profile (Fig. |^, solid line) fits 
the observed spectrum in the Q-branches of the 1/4 + 1/5 — 1/4 
band at 13.68 and 13.96 /im. Indeed, while the 1/5 band 
probes atmospheric levels between 2 and 5 mbar, the hot 
band sounds warmer, higher levels around 0.3 mbar, allow- 
ing us to determine the slope of the C2H2 profile between 
these two regions. 

Error bars on the C2H2 mixing ratio were estimated 
by taking into account instrumental noise and the uncer- 
tainty in the relative strengths of the C2H2 lines. The 
resulting mixing ratios are q — (8.9^.q[q) x 10~ 7 at 0.3 
mbar and q = (1-lio 1) x 10 -7 a t 4 mbar, giving a slope 
— d\nq/dhiP = 0.8 ± 0.1. The error on the temperature 
profile (±2K) introduces an additional uncertainty on the 
C2H2 mixing ratios of about 20%. This error, however, 
essentially equally applies to all pressure levels, and thus 
leaves the retrieved C2H2 profile slope mostly unaffected. 



3.2. Acetylene 

The ISO-SWS observations in the vicinity of the C2H2 
emission at 13.7 /im were compared with synthetic spectra 



3.3. Ethane 

Similarly to C2H2, we compare in Fig. || the ISO-SWS 
spectrum in the C2H6 1/9 band with synthetic spectra cal- 
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ciliated with three distinct C2H6 vertical profiles (Fig. |J). 
Each of the profiles was designed to reproduce the ob- 
served emission in the R Qo multiplet at 12.16 /mi, which 
probes the 1-mbar pressure level. The rest of the vc, band 
consists of weaker multiplets, which sound deeper levels 
extending from 1 to 10 mbar. This combination of strong 
and weak multiplets makes this band sensitive to the C2H6 
vertical distribution. In addition, the pseudo-continuum 
level in between the C2H6 emissions is also sensitive to 
the C2Hg abundance in the lower stratosphere. 

Our best-fit model (Fig. 0, solid line) has a slope 
— dlnq/d\nP = 0.6, but the steep-slope model (Fig. ||, 
dotted line) is also marginally compatible with the ob- 
servations. This results in a relatively large uncertainty 
011 the slope determination: q = (1.0 ± 0.2) x 10~ 5 at 
1 mbar, and q — (2.6+o'g) x 10 -6 at 10 mbar, giving 
— d\nq/d\nP — 0.6 ± 0.2. An additional error of 25% on 
q comes from the uncertainty on the temperature. Again, 
it should not affect the retrieved slope. 

3.4- Methylacetylene and diacetylene 

The ISO-SWS spectrum exhibits a broad emission near 
15.8 fim, which coincides with the ug band of CH3C2H 
(Fig. ||). Since the CH3C2H lines are optically thin, no in- 
formation on the vertical profile can be derived. Using a 
ver tical p rofile similar to that calculated by Gladstone et 
al. Q1996D , we found a column density of (1.5 ± 0.4) x 10 15 
molecule cm -2 . The synthetic spectrum exhibits small- 
scale structures which are not seen in the observations. 
This mismatch is propably due to an imperfect data re- 
duction. Indeed, the respective frequencies of the fringes 
and of the CH3C2H features lie close to each other. It is 
therefore very difficult to fully remove the former without 
altering the latter. On the contrary, the broad emission 
is a low frequency signal and is therefore left unaffected 
by the frequency filtering. This explanation is admittedly 
not entirely satisfactory, but given the wavelength match 
of the emission with the vg mode of CH3C2H, and in the 
absence of any other plausible candidates, we regard the 
detection of methylacetylene as unambiguous. 

The band of C4H2 is not detected at 15.92 /im. We 
inferred an upper limit of the C4H2 column density of 
7 x 10 13 molecule cm -2 , using the Gladstone et al. (1996) 
profile. 



4. Discussion 



Bezard et al. ( |1995[ ) first showed from 13.4-/im high- 
resolution spectroscopy that the C2H2 mixing ratio in- 
creases with height in the stratosphere, and that most of 
the acetylene is concentrated above the ~0.5-mbar level. 
Their distribution has a mixing ratio of about 7 x 10~ 7 
at 0.3 mbar, in reasonable agreem ent w ith our results. 
More recently, Betremieux & Yelle ( 199E ), using UV ob- 
servations, found an average C2H2 mixing ratio in the 20- 
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Fig. 5. Comparison between the ISO-SWS spectrum 
(solid line) and two synthetic spectra; with (dash-dotted 
line) and without (dotted line) CH3C2H and C4H2 opac- 
ities. 



60 mbar range of 1.5 x 10~ 8 , consistent with our value 
of 1.9 x 10~ 8 , extrapolated to this pressure range. Us- 
ing height-dependent mixing ratio profiles to analyse high- 
resolution infrared observations, Sada et al. (1998) derived 
mixing ratios of 3.9l* 3 x 10~ 6 for C2H6 at 5 mbar and 
2.3 ± 0.5 x 10" 8 for C 2 H 2 at 8 mbar. While the former 
value exactly agrees with our results, the latter is almost 
3 times less than that extrapolated downwards from our 
C2H2 profile. Also from infrared observations, Noll et al. 
( 1986[ ), using a slope of — dlnq/dlnP = 0.3, found a C2H6 
mixing ratio of 7.5 x 10 -6 at 1 mbar, which compares rea- 
sonably well with our measurement at the same pressure 
level. 

We also compared our retrieved mixing ratios with 
those calculated in the photochemical model of Gladstone 
et al. ( |1996[ ). For C2H6, our results are in excellent agree- 
ment with their model both at 1 and 10 mbar (Fig. ||). 
For C2H2, while the agreement is good at 4 mbar, their 
mixing ratio at 0.3 mbar is higher than ours by a factor of 
4. Our C2H2 slope is then significantly lower than theirs 
(- din qj din P = 0.8 ± 0.1 vs. -d\nq/d\nP = 1.2). Note 
that our derived C2H2 slope is still steeper than that of 
C2H 6 (— d\nq/d\nP = 0.6 ± 0.2), as expected. Indeed, 
C2H6, being less subject to photolysis, has a longer life- 
time in the jovian stratosphere than C2H2 (Gladstone et 
al. [1996| , their Fig. 6). 

Hydrocarbons are formed from the photolysis of CH4 
around the homopause. Small-scale turbulence, parame- 
terized in a photochemical model by the eddy diffusion 
coefficient (K), transports them downwards in the strato- 
sphere. This process is approximately modelled for long- 
lived products by the equation K(z)n(z)dq/dz — P(z), 
where n(z) is the number density at altitude z and P(z) 
the vertically integrated net production rate above alti- 
tude z. A first hypothesis would attribute the difference 
between the observed and calculated C2H2 abundances to 
an overestimation of the C2H2 production rate P{z) in 
the Gladstone et al. model. Our two abundance measure- 
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ments at 4 and 0.3 mbar allow us to calculate the mean 
dq/dz over this pressure range. Comparing with the value 
of dq/dz predicted by Gladstone et al. at 1 mbar, we found 
that P(z) should be decreased by a factor of ~3. 

A second explanation would be that Gladstone et al. 
underestimated the eddy diffusion coefficient K by a fac- 
tor of ~3 around 1 mbar. This underestimation at 1 mbar 
should also apply to the level and eddy diffusion coeffi- 
cient (Kh) at the homopause. However, it is difficult to 
directly quantify the changes on the homopause parame- 
ters, because of the strong coupling between production 
rates and homopause level. Simply, we note that our anal- 
ysis could imply an increase in Kh- It is consistent with 
Drossart et al. ( |1999| ), who, analysing CH4 fluorescence, 
found Kh = (7± 1) x 10 6 cm -2 s _1 , higher than the value 
of 1.4xl0 6 assumed in the Gladstone et al. model. In this 
case, the agreement on the C2H6 slope would also imply 
that the C2H 6 prod uction rate has been underestimated 
by Gladstone et al. (1996). In fact, the most direct con- 
clusion of our measurements is that Gladstone et al. have 
overestimated the C2H 2 /C2H 6 production rate ratio. 

The ISO-SWS spectrum enables the first detection of 
CH3C2H in the equatorial region of Jupiter. We retrieved 
a column density of (1.5±0.4) x 10 15 molecule cm -2 . Kim 
et al. (1985) had previously detected CH3C2H in the north 
polar auroral zone of Jupiter, and had retrieved a column 
density of (2.8± 2 ;i) x 1Ql6 molecule cm 2 . At least a large 
part of the difference is explained by different modelling 



assumptions. Kim et al. (1985) assumed a uniform vertical 



distribution throughout the stratosphere and used a tem- 
perature profile for the auroral regio n which is now known 
to be incorrect (Drossart et al. 1993 ). For C4H2, we found 
only an upper limit of 7 x 10 13 molecule cm -2 , 65 times 
lower than the stratospheric column density predicted by 
Gladstone et al. (4.5 xlO 15 molecule cm -2 ). This is con- 
sistent with their overestimate of C2H2 production, since 
the production of C4H2 is essentially quadratically depen- 
dent on the abundance of C2H2. As our CH3C2H column 
density is 3.5 times larger than calculated by Gladstone et 
al., we derived a CH3C2H/C4H2 ratio larger than 20, while 
they found a ratio of ~ 2. This discrepancy is all the more 
remarkable as, in the case of Saturn, wh ere bo th CH3C2H 
and C4H2 are detected (de Graauwet al. 1997 ), the photo- 
chemical model of Moses et al. ( 2000 ) gives, as observed, a 
CH3C2H/C4H2 ratio of about 10. This stresses that the C3 
and C4 chemistry in Jupiter should be reassessed with a 
complete photochemical model of the jovian stratosphere. 
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